Probing in-medium spin-orbit potential with intermediate-energy heavy-ion collisions 
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Incorporating for the first time both the spin and isospin degrees of freedom explicitly in transport 
model simulations of heavy-ion collisions at intermediate energies, we observe that a strong local 
spin polarization appears in the early stage of the reaction. Most interestingly, it is found that 
the nucleon spin up-down differential transverse flow is a sensitive probe of the spin-orbit potential 
independent of its tensor terms, providing a novel approach to probe both the density and isospin 
dependence of the in-medium spin-orbit coupling that is important for understanding the structure 
of rare isotopes and synthesis of superheavy elements. 
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The spin-orbit coupling is common to the motion of 
many objects in nature from quarks, nucleons, and elec- 
trons to planets and stars. In nuclear physics, it has been 
very well known for a long time that the spin-orbit cou- 
pling is crucial for understanding the structure of finite 
nuclei, such as the magic numbers [2, 0. However, many 
interesting questions regarding the in-medium spin-orbit 
coupling, especially its density and isospin dependence, 
remain unresolved although in free space it has been 
well determined from nucleon-nucleon (NN) scattering 
data [jj. Indeed, several phenomena observed or pre- 
dicted in studying properties of nuclear structures have 
been providing us some useful information about the in- 
medium spin-orbit interaction. For instance, the kink of 
the charge radii of lead isotopes can only be explained by 
introducing a weakly isospin-dependent spin-orbit cou- 
pling [1, HJ. Moreover, strong experimental evidences of 
a decreasing spin-orbit coupling strength with increas- 
ing neutron excess were reported 0, 0j- Furthermore, 
new experiments are currently being carried out at sev- 
eral laboratories to explore the density and isospin de- 
pendence of the spin-orbit coupling by comparing energy 
splittings of certain orbits in the so-called "bubble" nu- 
clei [1, Q with those in normal nuclei • The knowl- 
edge of the in-medium spin-orbit interaction is useful for 
understanding properties of drip-line nuclei [ll|, the as- 
trophysical r-process [l2j], and the location of stability 
island for superheavy elements n the other 

hand, it is also well known that the tensor force may af- 
fect the single-particle energy or even modify the magic 
numbers [15| . and may change the evolution of energy 
gaps between certain orbits with isospin asymmetry of 
finite nuclei [lij]. Moreover, calculations going beyond 
mean-field approximations indicate that the tensor force 
is also important for understanding short-range NN cor- 
relations [171 UM an d the density dependence of nuclear 
symmetry energy 19|, 2(j. Within Hartree-Fock calcula- 



tions, the spin-orbit interaction leads to a density gradi- 
ent term in the spin-orbit potential, while both the finite- 
range and tensor interactions contribute to the spin-orbit 
potential through the spin-current density terms, i.e., the 



tensor terms 21-23[. This mixing makes it very chal- 
lenging to disentangle respective contributions from the 
spin-orbit interaction and the tensor force to phenomena 
found in nuclear structure studies. It will be extremely 
useful if one could find ways to study them separately 
and determine their respective strengthes. In this Letter, 
within a newly developed spin-isospin dependent trans- 
port model, we show that the nucleon spin up-down dif- 
ferential transverse flow in heavy-ion collisions at inter- 
mediate energies is a sensitive probe of the in-medium 
spin-orbit potential independent of its tensor terms. 

Considering only time-even components of nuclear po- 
tentials, the single-particle Hamiltonian in nuclear mat- 



ter can be written as [24 1 



2m 



V 2 + U q + W q • (-»V x a). (q = n,p) (1) 



In the above, the second and third terms are the cen- 
tral and spin-orbit potential, respectively. The time-odd 
components are neglected as including them would gen- 
erally enhance dissipations 25j 26 1 but would not change 
the conclusion of the present work. Taking into account 
the density 27 [ and isospin dependence [1] as well as the 
tensor contribution [2l| , the form factor of the spin-orbit 
potential can be generally written as 



(aVp q + bVp q ,) + aJ q +(3J q ,, (q ^ q') (2) 



where p n (p) is the neutron (proton) density. W$(p) = 
Wo(p/ po) 1 represents the density-dependence of the spin- 
orbit coupling with p — p n +p p and Wo being the strength 
at saturation density po. Different values of 7, a, and 
b can be used to mimic its density and isospin depen- 
dence. We notice that with 7 = 0, a = 2 and 6=1 
Eq. (JSJ) reduces to the spin-orbit potential of the stan- 
dard Skyrme-Hartree-Fock model. While equal values 
for a and 6, and a nonzero value for 7 were predicted 
for the effective spin-orbit potential within a relativis- 
tic mean-field model [j| . The spin-current density J is 
the antisymmetric vector part of the spin-current tensor 
J^iv [HI- Neglecting the density dependence of Wq, the 
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spin-orbit coupling constant Wq ranges from about 100 
MeVfm 5 to 200 MeVfm 5 considering the uncertainty of 



the tensor terms [2l|, while the values of a and /3 rep- 
resenting the strength of the tensor terms are still quite 
uncertain. 

It is interesting to note that effects of the spin-orbit 
coupling were found earlier clearly in heavy-ion collisions 
at both low energies and relativistic energies. For in- 
stance, within the time-dependent Hartree-Fock (TDHF) 
calculations the nucleon spin-orbit potential was found 
to affect the fusion threshold energy [28| and leads to a 
local spin polarization 2(| 29[. In non-central relativis- 



tic heavy-ion collisions at 200 GeV/nucleon the partonic 
spin-orbit coupling was found to lead to a global quark 
spin polarization [3QJ . As the beam energy increases 
from low to intermediate energies around a few hundred 
MeV/nucleon, while effects of the spin-orbit coupling are 
expected to grow due to the higher angular momentum 
available, effects of NN scatterings which are mostly Pauli 
blocked at low energies also become increasingly impor- 
tant. Depending on how the spin of each particle changes 
during NN scatterings, the spin polarization might be 
strongly affected by the latter in heavy-ion collisions at 
intermediate energies. 

To model the spin-isospin dynamics in heavy-ion colli- 
sions at intermediate energies, we incorporate explicitly 
the spin degree of freedom and the spin-orbit potential 
in a previously developed isospin-dependent Boltzmann- 
Uchling-Uhlenbeck (BUU) transport model, see, e.g., 



Refs. [31|,|32J. The new model is dubbed SI-BUU12. To 
our best knowledge, the spin degree of freedom was never 
considered before in any of the existing transport models 
for heavy-ion reactions at intermediate energies since the 
emphasis of the community has been on extracting infor- 
mation about the Equation of State (EOS) of symmetric 
nuclear matter, density dependence of nuclear symme- 
try energy and in-medium NN scattering cross sections 
using spin-averaged experimental observables. In the SI- 
BUU12 model, the nucleon density p, spin density s, and 
the spin-current density J are all calculated by using the 



test particle method [331 . |34| . The equations of motion 



with the spin-orbit potential are 



dr 

db 
dp 

H 

da 
~dl 



P 



+ a x W a 



-VU q - V 



W q -(px a) 



— = 2(W q Xp) X <7. 



(3) 
(4) 
(5) 



In the center-of-mass (c.m.s) frame, two Fermi spheres 
of projectile/target nucleons are Lorentz boosted in the 
±z direction. During the collision process, the form fac- 
tor W points mainly along the impact parameter of the 
reaction, i.e., the x axis. Due to the spin-orbit poten- 
tial, the nucleon spin a tends to be parallel to p x W 
in order to lower the energy of the system. Namely, 



the nucleon spin is more likely to be polarized along the 
y axis perpendicular to the reaction plane (x-o-z). We 
will refer in the following text a nucleon with its spin 
in the +y (-y) direction as a spin-up (spin-down) nu- 
cleon. While our qualitative conclusions in this work 
are insensitive to the specific choice for U q , numerical 
results presented in the following are obtained by using 
a momentum-independent potential corresponding to a 
soft EOS of incompressibility Kq — 230 MeV, a symme- 
try energy E sym — 30 MeV and its density slope L = 60 
MeV at saturation density op = 0.16 fm -3 similar to the 
MSL0 parameterization 



Nucleon-nucleus scattering experiments have shown 
that nucleons may flip their spins after NN scatterings 
due to spin-related nuclear interactions, see, e.g., Ref. [36j 
for a review. Although not well determined yet, the spin- 
flip probability for in-medium NN scatterings is known to 
be appreciable, depending on the collision energy and the 
momentum transfer [37] • In the present work we will test 
different options of setting spins, including randomizing, 
flipping, or keeping spins unchanged after each NN scat- 
tering to study effects of different spin-flip probabilities 
on final observables. In addition, a spin- and isospin- 
dependent Pauli blocking is introduced in the SI-BUU12 
model. 
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FIG. 1: (color online) Contours of nucleon reduced density 
p/po (first row), y component of spin density s y (second row), 
and x component of the density gradient (Vp) x (third row) 
and spin-current density J x (fourth row) in the reaction plane 
at different stages in Au+Au collisions at a beam energy of 
400 MeV/nucleon and an impact parameter of b = 8 fm. Note 
that the x coordinates of the projectile (target) is negative 
(positive). 

As an illustration of spin-orbit coupling effects, it is 
interesting to first examine for a typical reaction at inter- 
mediate energies the time evolution of the y component 
of the spin density s y , the x component of the density 
gradient (Vp) x and the spin-current density J x in com- 
parison with the nucleon density contour in the reaction 
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plane. Shown in Fig. Q] are these quantities for Au+Au 
collisions at a beam energy of 400 MeV/nucleon and an 
impact parameter of b = 8 fm. For this example, we set 
Wo = 130 MeVfm 5 , 7 = 0, a = 2, b = 1, a = /3 = 0, 
and the spins of the colliding nucleons are randomized 
after each scattering. We will only change one or two 
of the above parameters while keeping all the rest the 
same in the following. Initially, we put the projectile and 
target nuclei without spin polarization far away so that 
a stable local spin polarization due to the spin-orbit po- 
tential can be generated in the two nuclei before they 
collide in the c.m.s frame. The initial distance is set to 
be long enough so that its further increment will not af- 
fect much the spin polarization of the two nuclei when 
they get in touch around 60 fm/c. Indeed, a local spin 
polarization is induced in the early stage of the reaction 
as expected. This "spin twist" phenomenon was first ob- 
served in TDHF calculations [26|, [29j . We note that the 
induced spin polarization observed here is stronger than 
that at lower energies due to the higher nucleon momen- 
tum at intermediate energies. It is seen that participant 
nucleons are more likely to spin-down while the specta- 
tor nucleons are more likely to spin-up during the early 
stage of the reaction. In the later stage, however, the spin 
polarization becomes weaker because of NN scatterings 
and other spin mixing effects, especially for participant 
nucleons at higher densities. At the end of the reaction, 
except for some large spectator-like fragments, the local 
spin polarization vanishes for most nucleons. 

As shown in the equations of motion, the spin-orbit 
coupling also affects the nucleon momentum and spatial 
distributions. In particular, spin-up and spin-down nu- 
cleons with the same momentum experienced opposite 
spin-orbit potentials during the whole collision process. 
Nucleon transverse collective flow, measured by using the 
average transverse momentum < p x (y) > in the reaction 
plane versus rapidity y, is one of the best known observ- 
able for revealing effects of density gradients in nuclear 
reactions (H)]. What can we learn about the spin-orbit 
potential from the nucleon transverse flow? To answer 
this question, first we compare in the left panel of Fig. [2] 
the transverse flows of spin-up and spin-down nucleons. 
It is seen that the transverse flow of spin-up nucleons is 
larger than that of spin-down ones. This can be easily 
understood by looking at the x-component of the den- 
sity gradient shown in the third row of Fig. [TJ By ex- 
amining the time evolution, we found that the effects of 
the spin-orbit potential on the transverse flow during the 
first 20 fm/c of the collision are mostly washed out due 
to violent interactions. The spin-dependent transverse 
flow is mainly determined by the dynamics afterwards. 
As the projectile (target) is still moving in the +z (- 
z) direction, the participant nucleons from the projectile 
(target) with negative (positive) (Vp)^ give a more repul- 
sive/attractive spin-orbit potential [W ■ (p x a) > / < 0] 
for spin- up/down nucleons. This leads to a larger trans- 
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FIG. 2: (color online) Transverse flow of spin-up nucleons and 
spin-down nucleons (a) and spin up-down differential flows 
using different Wo (b) for the same reaction as in Fig.l. 



verse flow for spin-up nucleons than spin-down ones. 

To extract more accurately information about the spin- 
orbit potential without much hinderance of central po- 
tentials, we investigate next the spin up-down differential 
transverse flow 
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where N(y) is the number of nucleons with rapidity y, 
and <7i is 1(— 1) for spin-up (spin-down) nucleons. Simi- 
lar to the neutron-proton differential transverse flow for 
probing the symmetry potential 39], the spin up-down 
differential transverse flow maximizes effects of the oppo- 
site spin-orbit potentials for spin-up and spin-down nu- 
cleons while canceling out largely spin-independent con- 
tributions. Indeed, the spin up-down transverse flow is 
a sensitive probe of the spin-orbit coupling strength Wq. 
As an example, shown in the right panel of Fig. [2] is a 
comparison of the spin up-down differential transverse 
flows obtained using the two limiting values of Wq used 
in the literature. As expected, a larger value of Wq gives 
a stronger spin-orbit potential and thus a larger spin up- 
down differential transverse flow. 

The density dependence of the spin-orbit coupling, 



which was tested earlier, see, e.g., Ref. [27|, is still al- 
most completely unknown, and this has recently moti- 
vated more new experiments [loj ]. To investigate effects 
of the density dependence of the spin-orbit coupling on 
the spin up-down differential flow, shown in the panel 
(a) of Fig. [3] are the results obtained by varying only the 
7 parameter. It is seen that the spin up-down differen- 
tial flow is larger for a weaker density dependence of the 
spin-orbit coupling if the strength at saturation density 
is fixed. 

The isospin dependence of the spin-orbit coupling is 
another interesting issue especially relevant for under- 
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FIG. 3: (color online) Spin up-down differential transverse 
flows using different density dependence of the spin-orbit cou- 
pling coefficients (a) and from pure like-nucleon (b) and pure 
unlike-nucleon coupling (c) for the same reaction as in Fig.l. 
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FIG. 4: (color online) Comparison of the spin up-down differ- 
ential transverse flows using different coefficients of the tensor 
terms (a) and different treatments of spin after scatterings (b). 



standing the structure of rare isotopes and the synthesis 
of superheavy elements. We notice that available studies 
on nuclear structure seem to favor a weak isospin depen- 
dence of the spin-orbit coupling 0,0, EH- To evaluate po- 
tential applications of our approach in further constrain- 
ing the isospin dependence of the spin-orbit potential, 
we next compare the spin up-down differential flows for 
neutrons and protons using the pure like-nucleon cou- 
pling (a = 3 and 6 = 0) and pure unlike-nucleon cou- 
pling (a — and b — 3) in panel (b) and panel (c) of 
Fig. [3l respectively. As the system considered is globally 
neutron-rich and Vp n is generally larger than V/3 P , the 
pure like(unlike)-nucleon coupling leads to an apprecia- 
bly larger (smaller) spin up-down differential flow for neu- 
trons than for protons. Moreover, the unlike-nucleon cou- 
pling generally reduces significantly the overall strength 
of the spin-orbit potential and thus the spin up-down dif- 
ferential flow. Of course, more neutron-rich systems will 
be better for probing the isospin-dependence of the spin- 
orbit coupling using the double differential flow between 
spin up-down neutrons and protons. 

So far we haven't included in the spin-orbit potential 
contributions from the tensor terms yet. With spin po- 
larization, as a is parallel to px Vp, J ~ px a is generally 
opposite to Vp. Thus, J terms with positive (negative) 
coefficients contribute negatively (positively) to W, as 
mentioned in Ref. [29j . Although great efforts have been 
made in determining the values of a and j3 from studying 



the structure of finite nuclei [2lM23| . they are still quite 
uncertain. By setting a = /3 = 100 or —100 MeVfm 5 , 
which are considered reasonable in magnitude based on 
the early studies in the literature, we compare the spin 
up-down differential transverse flow in the left panel of 
Fig. 01 It is seen clearly that the spin up-down differen- 
tial transverse flow is rather insensitive to the J terms. 
This is understandable as the Vp terms related to the 



density distribution remain strong while the spin-current 
density vanishes after the first 20 fm/c of the collision as 
shown in Fig. [T] 

Finally, what are the effects of the possible spin flip 
in NN scatterings on the spin up-down differential flow? 
We answer this question by using the results shown in 
the right panel of Fig. [4] Due to the lack of knowl- 
edge about the energy and isospin dependence of the in- 
medium spin-flip probability, we compare the following 
three choices for setting spins of nucleons after NN scat- 
terings: (1) flipped, (2) randomized, or (3) unchanged, 
effectively from a large to small spin-flip probability. It 
is seen that the spin up-down differential transverse flow 
decreases with the increasing spin-flip probability as one 
expects. Moreover, it is very encouraging to see that the 
spin up-down differential flow is still considerable even if 
a 100% spin-flip probability is assumed. 

In summary, the spin degree of freedom and the spin- 
orbit potential are incorporated for the first time in 
an isospin-dependent transport model providing a useful 
new tool for investigating the spin-isospin dynamics of 
heavy-ion collisions at intermediate energies, such as the 
development of local spin polarization in these reactions. 
The nucleon spin up-down differential transverse flow is 
shown to be a sensitive probe of the in-medium spin-orbit 
potential independent of its tensor terms. Comparisons 
with future experiments will allow us to determine the 
density and isospin dependence of the in-medium spin- 
orbit coupling that has significant ramifications in both 
nuclear physics and astrophysics. 
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